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The main properties of realistic models for manganites are studied using analytic mean-field 
approximations and computational numerical methods, focusing on the two-orbital model with elec- 
trons interacting through Jahn- Teller (JT) phonons and/or Coulombic repulsions. Analyzing the 
model including both interactions by the combination of the mean-field approximation and the 
exact diagonalization method, it is argued that the spin-charge-orbital structure in the insulating 
phase of the purely JT-phononic model is not qualitatively changed by the inclusion of the Coulomb 
interactions. As an important application of the present mean-field approximation, the C-E-type 
AFM state, the charge-stacked structure along the z-ajds, and the {Zx^ — r^)/{Sy^ — r^)-type orbital 
ordering are successfully reproduced based on the JT-phononic model for the half-doped manganite, 
in agreement with recent Monte Carlo simulation studies by Yunoki et al. Topological arguments 
and the relevance of the Heisenberg exchange among localized t2g spins explains why the inclusion 
of the nearest-neighbor Coulomb interaction does not destroy the charge stacking structure. It is 
also verified that the phase-separation tendency is observed both in purely JT-phononic and purely 
Coulombic models in the vicinity of the hole undoped region, as long as realistic hopping matrices 
are used. This highlights the qualitative similarities of both approaches, and the relevance of mixed- 
phase tendencies in the context of manganites. In addition, the rich and complex phase diagram of 
the two-orbital Coulombic model in one dimension is presented. Our results provide robust evidence 
that Coulombic and JT-phononic approaches to manganites are not qualitatively different ways to 
carry out theoretical calculations, but they share a variety of common features. 



Andre Luiz Malvezzi 
Departamento de Fisica - Faculdade de Ciencias - Universidade Estadual Paulista 
Caixa Postal 473, 17.033-360, Bauru, SP, Brazil 



PACS numbers: 71.10.-w, 75.10.-b, 75.30.Kz 



The recent discovery of the colossal magnetoresistance 
(CMR) phenomenon in manganese oxides has triggered 
a huge exnpi-imcntal and theoretical effort to understand 
its origin.EfH The large changes in resistivity observed in 
experiments usually involve both a low-temperature or 
high magnetic field ferromagnetic (FM) metallic phase, 
and a high-temperature or low magnetic field insulating 
phase. For this reason, to address the CMR effect, it ap- 
pears unavoidable to have a proper understanding of the 
phases competing with the metallic FM regime, which as 
a first approximation can be rationalized based on the 
standard double exchange (DE) mechanism.&a In fact, 
the phase diagram of manganites has revealed a very 
complex structure, clearly showing that metallic ferro- 
magnetism is just one of the several spin, orbital, and 
charge jaxrangements that are possible in the manganese 
oxides. B'EI 



I. INTRODUCTION 



culations involved only one orbital, and the many-body 
techniques used in its analysis were mostly restricted 
to crude mean-field approximations (MFA). Quite re- 
cently, the first computational studies of the one-orbital 
model £pr manganites have been presented by Yunoki 
et aLLrEI While the expected FM-phase emerged clearly 
from such analysis, several novel features were identified, 
notably phase-separation (PS) tendencies close to n—1 
(where n is the Cg electron number density per site) be- 
tween an antiferromagnetic (AFM) insulating phase and 
a metallic FM state. Several calculations have confirmed 
these tendencies toward mixed-phase characteristics, us- 
ing a variety of techniques. t2l Moreover, the computa- 




The theoretical study of manganites started decades 
ago when the DE ideas to explain the FM-phase were pro- 
posed. Unfortunately, the model used in those early cal- 



spin-FM phases that differ in their orbital arrangement, 
i.e., staggered vs. uniform, were again identified using 
one-dimensional (ID) and two-dimensional (2D) clusters. 
This illustrates the relevance of the orbital degree of free- 
dom in the analysis of manganites, and confirms the im- 
portance of mixed-phase characteristics in this context. 
The theoretical, mostly computational, calculations have 
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been recently summarized by Moreo et al. ,0 where it has 
been argued that a state with clusters of one phase em- 
bedded into another should have a large resistivity and a 
large compressibility. Other predictions in mixed-phase 
regimes involve the ^presence of a robust pseudogap in 
the density of states Jlj similar to the results obtained us- 
ing angle-resolved nhptoemission experiments applied to. 
bilayer manganitesJlj A large number of experimentsE£l 
have reported Jie presence of inhomogeneities in real 
manganites,EjilZl results compatible with those obtained 
using computational studies. More recently, the influence 
of disorder on metal-insulator transitions of manganite 
models that would be of fi*st-order character without 
disorder has been proposedlla as the origin of the large 
cluster coexiS|teiice reported in experiments for several 
compounds.t^llZl It is clear that intrinsic tendencies to- 
ward inhomogeneous states is at the heart of manganite 
physics, and simple scenarios involving small polarons are 
not sufficient to understand the mixed-phase tendencies 
observed experimentally. 

Although the presence of mixed-phase tendencies com- 
peting with the FM-phase is by now well-established and 
a considerable progress has been achieved in the theo- 
retical study of models for manganites, many issues still 
remain to be investigated in this context. One of them is 
related with the notorious charge ordering phenomenon, 
characteristic of the narrow-band manganites. Especially 
at n=0.5, the so-called Ci?-type AFM phase has been 
widely observed in materials such as Lao.sCao.sMnOa 
and Ndo.sSrg.sMnOs. In spite of the considerable impor- 
tance of this phase in view of its relevance for the OMR ef- 
fect reported at high hole density,B it is only recently that 
the Ci?-phase has been theoretically undeps.taDd using 
the analytic MFA in the JT-phonoiui;pHiodeltj'L3 and nu- 
merical computational techniques Ed'c3 The importance 
of the zigzag chains of the Ci?-state, along which the 
spins are parallel, has been remarked in these recent 
efforts. In fact, it has been shown that even without 
electron-phonon or Coulomb interactions the Ci?-state 
structure is nevertheless stable, and in this limit it arises 
from a simple band-insulator picture. Previous-|Work, 
such as the pioneering results of Goodenough,E3 were 
based on the assumption of a checkerboard charge state 
upon which the spin and orbital arrangements were prop- 
erly calculated, but they did not consider the competition 
with other charge disordered phases in a fully unbiased 
calculation as carried out in Refs. [T9|-p^. 

Here a particular feature of the CO Ci?-phase struc- 
ture should be remarked. In the x-y plane, the Mn'^+ 
and Mn'*+ ions arrange in a checkerboard pattern, which 
naively seems to arise quite naturally from thp presence 
of a strong long-range Coulomb interactiono However, 
contrary to this naive expectation, the same CO pat- 
tern stacks in experiments along the z-axis without any 
change, although the t2g spin direction alternates from 
plane to plane. In other words, the Mn^+ and Mn''+ 
ions in the Ci?-type AFM structure do not form a three- 
dimensional (3D) Wigner-crystal structure, clearly sug- 



gesting that the origin of the CO phase in the half- 
doped manganite cannot be caused exclusively by a dom- 
inant long-range Coulomb interaction. Other interac- 
tions should be important as well. In this work, it is 
shown that the charge-stacked (CS) CO-state in the CE- 
type structure found in experiments can also be obtained 
in the simple MFA for the purely JT-phononic model. In 
addition, it is shown that this structure is not easily de- 
stroyed by the inclusion of the nearest-neighbor Coulomb 
interaction. The exchange Jaf between the i2g spins 
plays a key role in the stabilization of the CS structure. 
The details of the calculation are discussed here. The 
present semi-analytical results are in excellent agreement 
with recent Monte Carlo (MC) simulation^ and provide 
a simple formalism to rationalize these computational re- 
sults. 

On the other hand, it is possible to take another ap- 
proach to understand the CO state in manganites in- 
cluding the Ci?-type structure, by|-emphasizing the role 
of the on-site Coulomb interaction£j In this direction of 
study, several results are available in the literature. For 
instance, Ishihara et al. discussed the spin-charge-orbital 
ordering from the limit of infinite strong electronic cor- 
relation, nacnely, under the approximation of no double 
occupancyES Maezono et al. provided information on the 
overall features of the phase diagram for manganites by 
using a mean-field calculation.c3 These two mean-field 
approaches do not reproduce the Ci?-state of half-doped 
manganites. Mizokawa and Fujimori discussed the sta- 
bilization of the Ci?- type stmcture by using the model 
Hartree-Fock approximation.LZl However, to the best of 
our knowledge, the origin of the CS structure in the CE- 
type state has not been clarified based on the Coulombic 
scenario, although several properties of manganites have 
been reproduced. In this paper, a possible way to un- 
derstand the Ci?-type state with the CS structure in the 
purely Coulombic model even without the long-range re- 
pulsion is briefly discussed. 

The present paper contains information about other 
subjects as well. In spite of the importance of additional 
interactions besides the long-range Coulombic one to sta- 
bihze the proper charge-stacked n=0.5 state, as discussed 
before, it is clear that at least the on-site Coulomb inter- 
actions are very strong and their influence should be con- 
sidered in realistic calculations. Most of the previous an- 
alytic and computational works that reported the PS ten- 
dency and the CO phase have used electrons interacting 
among themselves indirectly through their coupling with 
localized i2g spins or with JT phonons. The JT-phononic 
model has been extensively studied by the present au- 
thors since the explicit inclusion of on-site Coulomb inter- 
actions diminish dramatically the feasibility of the com- 
putational work. In addition, it has been found that 
the JT-phononic model explains quite well several exper- 
imental results, even if the Coulomb interaction is not in- 
cluded explicitly. Moreover, it should be noted that the 
energy gain due to the static JT distortion is maximized 
when only one Cg electron is located at the JT center. 
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Confirming this expectation, the double-occupancy of a 
given orbital has been found to be negligible in previ- 
ous investigations at intermediate and large values of the 
electron-phonon coupling, and in this situation, adding 
an on-site repulsion would not alter the physics of the 
state under investigation. Nevertheless, although the 
statements above are believed to be qualitatively correct, 
it should be checked explicitly in a more realistic model 
that indeed the physics found in the extreme case of only 
phononic interactions survives when both JT-phononic 
and Coulombic interactions are considered. Thus, an- 
other purpose of this paper is to clarify the validity of 
the JT-phononic only model by showing that the inclu- 
sion of the Coulomb interaction does not bring qualita- 
tive changes in the conclusions obtained from the purely 
phononic model. For this purpose, an analytic-numeric 
combination is employed, namely, the MFA including 
Coulombic terms is developed in detail and the results 
are compared with exact diagonalization (ED) data for 
a small cluster. In addition, the PS tendency is reex- 
amined briefly within the framework of MFA, and some 
comments on previous results are provided. It is clear 
that addressing whether purely JT phononic or purely 
Coulombic approaches lead or not to qualitatively differ- 
ent phase diagrams is an important area of investigation. 

As a special case of the general goal described in the 
previous paragraph, the issue of PS in the presence of 
Coulomb interactions will also be studied in this pa- 
per using computational techniques. In fact, it has not 
been analyzed whether models with two orbitals and 
Coulomb repulsions, i.e. without JT-phonons, lead to 
mixed-phase tendencies as pure phononic models do. In 
this paper it is reported that a study of a simple one- 
dimensional two-orbital model without phonons but with 
strong Coulomb interactions, produces PS tendencies 
similadjf as reported in previous investigations by our 
group.Ell Several features are in qualitative agreement 
with those observed using JT-phonons, illustrating the 
similarities between the purely phononic and the purely 
Coulombic approaches to manganites, and the clear rele- 
vance of mixed-phase characteristics for a proper descrip- 
tion of these compounds. It is concluded that simple 
models for manganites, once studied with robust compu- 
tational techniques, are able to reveal a complex phase 
diagram with several phases having characteristics quite 
similar to those observed experimentally. 

The organization of the paper is as follows. In Sec. II, 
a theoretical model for manganites is introduced, and 
several reduced Hamiltonians are discussed in prepara- 
tion for the subsequent investigations. Section III is de- 
voted to the MFA. The formulation for the MFA is pro- 
vided and a comparison with the ED result is shown 
to discuss the validity of the MFA. It is argued that 
the JT-phononic interaction plays a role more relevant 
than the Coulombic interaction as suggested from the 
viewpoint of the continuity of the orbital-ordered (00) 
phase. In Sec. IV, as an application of the present MFA, 
the CO/00 structure in half-doped manganites is stud- 



ied. Especially, the origin of the Ci?-type state with 
the CS structure is discussed in detail. In Sec. V, the 
PS tendency is studied in the framework of the MFA 
for the purely JT-phononic and Coulombic models. In 
Sec. VI, the two-orbital ID Hubbard model is analyzed 
using the density matrix renormalization group (DMRG) 
technique. Even in the purely Coulombic model, PS ten- 
dencies and CO phases are observed. Finally in Sec. VII, 
the main results of the paper are summarized. The main 
overall conclusion is that studies based upon purely JT- 
phononic or Coulombic approaches do not differ substan- 
tially in their results, and there is no fundamental differ- 
ence between them. For reasons to be described below, of 
the two approaches the best appears to be the phononic 
one due to its ability to select the proper orbital order- 
ing pattern uniquely. PS appears in purely phononic and 
purely Coulombic approaches as well. In Appendix A, 
the relation between the PS scenarios for cuprates and 
manganites is discussed. Throughout this paper, units 
such that h=kB=l are used. 



II. MODELS FOR MANGANITES 

In order to explain the existence of ferromagnetism 
at low temperatures in the mangaiiites, the DE mech- 
anism has been usually invoked.l3u In this framework, 
holes can improve their kinetic energy by polarizing the 
spins in their vicinity. If only one Cg orbital is_used, 
the FM Kondo or one-orbital model is obtained,L3 and 
several resiilis for this model are already available in the 
literatureUla Although the important PS tendencies have 
been identified in this context, it is clear that the one- 
orbital model is not sufficient to describe the rich physics 
of the manganites, where the two active orbitals play a 
key role and the CO state close to density 0.5 is crucial 
for the appearance of a huge MR effect at this density. 
Thus, in this paper, the two-orbital model is discussed to 
understand the manganite physics, by highlighting the 
complementary roles of the JT phonon and Coulomb in- 
teractions. 



A. Hamiltonian 

Let us consider doubly-degenerated Cg electrons, 
tightly coupled to localized t2g spins and local distor- 
tions of the MnOe octahedra. This situation is described 
by the Hamiltonian H composed of five terms 

H = iJkin + ^^Hund + ^^AFM + ^^cl-ph + ^^cl-ol- (1) 

The first term indicates the hopping motion of eg elec- 
trons, given by, 

ia77'(T 
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where diao- {diha) is the annihilation operator for an Cg- 
electron with spin a in the d^^-y^ (d^z^-r^) orbital at 
site i, a is the vector connecting nearest-neighbor sites, 
and t^^, is the nearest-neighbor hopping amplitude be- 
tween 7- and 7'-orbitals along the a-direction. The am- 
plitudes are evaluated from the overlap integral between 
manganese and oxygen ionspj given by 

iaa = -VSt^b = -A/3tL = 3tSb = t, (3) 

for the cc-direction, 

tL = Vstl^ = V3tL = 3i^, = (4) 
for the y-dircction, and 

t^b-4i/3,tL-t|b-C = 0, (5) 

for the z-direction. Note that t^^ (equal to t^^ by sym- 
metry) is taken as the energy scale t. 

The second term is the Hund coupling between the eg 
electron spin Si and the localized t2g spin Si, given by 

^^Hund = — >/h ^ Si • Sj, (6) 
i 

with Si = J^ja/s dl^a'^afsdijis, whcrc Jh(>0) is the Hund 
coupling, and cr=(ax, Uy, CTz) are the Pauli matrices. 
Note here that the t2g spins are assumed classical and 
normalized to |Si|=l. 

The third term accounts for the G-type AFM property 
of manganites in the fully hole doped limit, given by 

Hafm = Jaf ^ Si • Sj, (7) 

(i,j> 

where Jaf is the AFM coupling between nearest neigh- 
bor t2g spins. The value of Jaf will be small in units of 
t to make it compatible with experiments for the fully 
hole-doped CaMnOa compound.E] 

In the fourth term, the coupling, oiieg electrons to the 
lattice distortion is considered asEirE3 

^^ol-ph = gY^[Qli{dl^„diaa + dl^dita) 
ia 

+ Q2i{d\^„diha + d\^^diaa) + Qsiid^A^aa - dl^diba)] 

+ (1/2) ^[fcb.Q?i + kMQl + Qli)l (8) 

i 

where g is the coupling constant between Cg electrons and 
distortions of the MnOg octahedron, Qn is the breathing- 
mode distortion, Q2i and Qyi are, respectively, the JT 
distortions for the (x^ — y^)- and (3z^ — r^)-type modes, 
and fcbr (^jt) is the spring constant for the breathing 
(JT) mode distortions. Here the distortions are treated 
adiabatically, but the validity of this treatment should 
be discussed. In general, the adiabaticity is judged from 
the ratio rj— Tc/ti, where To and Tj are characteristic time 



scales for the motion of electrons and ions, respectively. 
If 77 is much less than unity, the adiabatic approximation 
holds, since the motion of ions is very slow compared 
to that of electrons. Due to the relations Tj^cj"^ and 
Tc^t~^, where to is the frequency of the lattice distor- 
tioTx—ri^Lu/t is obtained. In the manganite, t is 0.2-0.5 
eV,^ while lu is about 500-600cm-i.ij Thus, the adia- 
batic approximation is in principle valid in the study of 
the JT distortion of manganites. An important energy 
scale characteristic of i/ci-ph is the static JT energy, de- 
fined by 

EjT = 5V(2fc,TT), (9) 

which is naturally obtained by the scaling of Q^i = 
{g/kjT)q^i with /i=l,2, and 3, where g/kjT is the typical 
length scale for the JT distortion. By using Ejt and t, it 
is convenient to introduce the non-dimensional electron- 
phonon coupling constant A as 

A - ^2EiT/t. (10) 

The characteristic energy for the breathing-mode distor- 
tion is given by £'br=5^/(2fcbr)=£'jT//3 with (3=ki,^/kjT- 
The last term indicates the Coulomb interactions be- 
tween Cg electrons, expressed by 

Hc\-c\ = Pi7TPi7i + ^' ^ PiaaPiha' 

+ J '^Aaadlha'die.a'diha +Vj2PiPh (H) 
iaa' (ij) 

with pi^a= dj^^di^a- and pi=J2^„ Pi-ya, where U is the 
intra-orbital Coulomb interaction, U' is the inter-orbital 
Coulomb interaction, J is the inter-orbital exchange in- 
teraction, and V is the nearest-neighbor Coulomb inter- 
action. Note that the Hund-coupling term iJuund be- 
tween conduction and localized spins also arises from 
Coulombic effects, but in this paper, the "Coulomb in- 
teraction" refers to the direct electrostatic repulsion be- 
tween Eg electrons. Note also that the relation U —U'+2J 
holds in the localized ion system. Although the validity 
of this relation may not be guaranteed in the actual ma- 
terial, it is assumed to hold in this paper to reduce the 
parameter space in the calculation. 

Finally, a comment about the cooperative JT effect 
which is not included in the phononic models is discussed 
here. In principle, a Heisenberg-like coupling between 
nearest-neighbor phonons should also be added to ac- 
count for the fact that when a given octahedra is elon- 
gated iupOpc direction, its neighbors deform in an oppo- 
site way.E2l However, this cooperative effect has not been 
studied in detail using computational techniques and here 
the focus of the paper will be on results obtained with 
non-cooperative phonons. The rich phase diagram and 
good qualitative agreement with experiments justifies a- 
posteriori this extra approximation. In fact, in the un- 
doped limit, there are no essential differences between 
the results with and without the cooperative effect .E3 
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B. Reduced Hamiltonians 

The Hamiltonian Eq. (Q) is believed to define an ap- 
propriate model for manganites, but it is quite difficult 
to solve it exactly. In order to investigate further the 
properties of manganites based on H, some simplifica- 
tions and approximations are needed. A simplification 
without the loss of essential physics is to take the widely 
used limit Jh=oo. In such a limit, the Cg electron spin 
perfectly aligns along the t2g spin direction, reducing the 
number of degrees of freedom. Moreover, there is a clear 
advantage that both the intra-orbital on-site term and 
the exchange term can be neglected in _ffci-ci- Thus, the 
following simplified model is obtained: 

H°° = — ^ '5i4+ai^T.'cJ^Ci+a7' + JaF ^ Si • Sj 
ia77' (ij) 

+ -EjT ^[2(giini + (72iTxi + qsiTzi) + fiqli + qli + ql\\ 

i 

+ U'^ni^m, + V^nin-^, (12) 



where is the spinless e„ electron operator, given b; 



Ci^— cos(0i/2)di^|-|- sin(0i/2)e the angles 9i and 

(pi define the direction of the classical t2g spin at site 

i, ^T'i'y— Cj^Ci-y, '^\—^22t'y ^175 '^x\— Cj^Cifo-t-Cj^Cia, and Tz\ = 

cJ^Cia— Cj^Cifc. Here 5ij denotes the change of hopping 
amplitude due to the difference in angles between 
spins at sites i and j , given by 



5ij = cos(6ii/2)cos(6'j/2) 

+ sin(6'i/2)sin(6lj/2)e-*('^'-"*J 



(13) 



In principle, 9i and 0i should be optimized to provide 
the lowest-energy state. However, in the following ana- 
lytic approach, only spin configurations such that nearest 
neighbor i2g spins are either FM or AFM will be consid- 
ered. Namely, the spin canting state is excluded from the 
outset. The validity of this simplification is later checked 
by comparison with the numerical results obtained from 
a relaxation technique. Note that this simplification does 
not restrict us to only fully FM or G-type AFM states, 
but a wide variety of other states, such as the Ci?-type 
one, can also be studied. 

Based on the reduced Hamiltonian iJ °° , it is possible to 
consider two scenarios for manganites, as schematically 
shown in Fig. 1. Temporarily, V will not be taken into 
account to focus on the competition between £^jt and U' . 
One of them is the JT-phonon scenario in which U' is con- 
sidered as a correction into the purely JT-phononic model 
H^=H°°{U'=Qy=Q) taken as starting point. Another 
is the Coulomb scenario where the JT distortion is a cor- 
rection into the purely Coulombic Hubbard-like model 
Hq'=H°°{Ejt=0)- Of course, if the many-body analy- 
sis were very accurate then there should be no difference 
between the final conclusions obtained from both sce- 
narios when the realistic situation, £^jt>0 and U'>0, is 



approached. However, if this realistic situation can be 
continuously obtained from the limiting cases U'=0 or 
E]T—0, then it is enough to consider only the limiting 
models iJjf- or to grasp the essence of the mangan- 
ite physics. On the other hand, it may also occur that 
starting in one of the two extreme cases a singularity ex- 
ists preventing a smooth continuation into the realistic 
region of parameters. In this case, just one of the sce- 
narios would be valid. For investigations in this context, 
the simplified model H°° certainly provides a good stage 
to examine the roles of U' and i?jT- In the following 
section, this point will be discussed in detail. 
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FIG. 1. Schematic representation of the two approaches 
(Jahn- Teller and Coulombic) to the realistic situation based 
on the simplified model H°° . 

Another possible simplification could have been ob- 
tained by neglecting the electron-electron interaction in 
H but keeping the Hund coupling finite, leading to the 
following purely JT.-phononic model with active spin de- 
grees of freedomiOEj 



^^.IT — ^kin + ^^Hund + HafM + ^^cl-ph- 



(14) 



Note that this model is reduced to Hy^ when the limit 
Jh=oo is taken. To solve -ffjT, numerical methods such 
as MC techniques and-the relaxation method have been 
employed in the past.E3E3 However, here this model will 
not be discussed explicitly, but the simplified version 
Hj^ will be investigated in the next section. Qualita- 
tively, the negligible values of the probability of double 
occupancy at large A justifies the neglect of Hei-ei- In 
addition, the results of recent studies addressing the A- 
type AFM-phase of the hole-undoped limit using coop- 
erative JT-phonons have been shown to be quite similar 
to those found with pure Coulombic interactions,Lj the 
latter treated with the MFA. 

Nevertheless, in spite of the above discussed indica- 
tions that the JT and Coulomb formalisms lead to sim- 
ilar physics, it would be important to verify this belief 
by studying a multi-orbital model with only Coulombic 
terms, without the extra approximation of using mean- 
field techniques for its analysis. Of particular relevance 
is whether PS tendencies and charge ordering appear in 
this case, as they do in the JT-phononic model. This 
analysis is particularly important since, as explained be- 
fore, a mixture of phononic and Coulombic interactions 
is expected to be needed for a proper quantitative de- 
scription of manganites. For this purpose, yet another 
simplified model will be analyzed in this paper: 
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(15) 



Note that this model cannot be exactly reduced to H^, 
since the Hund coupling term between Cg electrons and 
t2g spins is not explicitly included. The reason for this ex- 
tra simplification is that the numerical complexity in the 
analysis of the model is drastically reduced by neglecting 
the localized t2g spins. In the FM phase, this is an excel- 
lent approximation but not necessarily for other magnetic 
arrangements. Nevertheless the authors believe that it 
is important to establish with accurate numerical tech- 
niques whether the PS tendencies are already present in 
this simplified two-orbital models with Coulomb interac- 
tions, even if not all degrees of freedom are incorporated 
from the outset. Adding the S=3/2 quantum localized 
spins to the problem would considerably increase the size 
of the Hilbert space of the model, making it intractable 
with current computational techniques. 



III. MEAN-FIELD APPROXIMATION FOR H°° 

Even a simplified model such as is still difficult to 
be solved exactly except for some special cases. Thus, in 
this section, the MFA is developed for H°° to attempt to 
grasp its essential physics. Note that even at the mean- 
field level, due care should be paid to the self-consistent 
treatment for the lift of the double degeneracy of eg elec- 
trons. 



A. Formulation 

First let us rewrite the electron-phonon term by ap- 
plying a simple standard decoupling such as q2iT3:i ~ 
{q2i)Txi + q2i{Txi) - {q2i){Txi), where the bracket means 
the average value using the mean-field Hamiltonian de- 
scribed below. By minimizing the phonon energy, the lo- 
cal distortion is determined in the MFA as qii=— (rii)//3, 
q2i=—{Txi), and q3i=—{Tzi). Thus, the electron-phonon 
term in the MFA is given by 

-ffcY^h = -2^[-Bbr(rii)ni + EjT{{Txi)Txi + (T,i)T,i)] 
i 

+ $^[i?br(ni)' + i?JT((T.i)2 + (T,i)2)]. (16) 
i 

Now let us turn our attention to the electron-electron 
interaction term. At a first glance, it appears enough to 
make a similar decoupling procedure for iJoi-ci- How- 
ever, such a decoupling cannot be uniquely carried out, 
since Hd-oi is invariant with respect to the choice of Cg- 
electron orbitals due to the local SU(2) symmetry in the 
orbital space. Thus, to obtain the 00 state, it is neces- 
sary to find the optimal orbital set by determining the 
relevant Cg-electron orbital self-consistently at each site. 
For this purpose, it is useful to express q2i and qsi in 
polar coordinates as 



q2i= qi sin , q3i= qi cos . 
In the MFA, qi and can be determined as 

qi - ViTxi)^ + (Tzi)^^ 

and 
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tan-i((T,i)/(T,i)). 



(17) 



(18) 



(19) 



By using the phase ^i, Cia and Cib are transformed into 
Cia and Cib as Qa = e'^'/^[cia cos(6/2) -I- Qb sin(^i/2)] and 
= e'^'/2[-CiaSin(^i/2) + Cib cos(Ci/2)]. Note that the 
phase factor is needed for the assurance of the single- 
valuedness of the basis function, leading to the molecular 
Aharonov-Bohm effect. E2l 

It should be noted that the phase determines the 
electron orbital set at each site. For instance, at ^i=27r/3, 
"a" and "b" denote the dy2_z2- and c?32.2_^2-orbital, re- 
spectively. In Table. I, the correspondence between 
and the local orbital is summarized for several important 
values of ^i. Note here that d33.2_^2 and d3y2_^2 nevp; 
appear as the local orbital set. In recent publications,t2l 
those were inadvertently treated as an orthogonal orbital 
set to reproduce the experimental results, but such a 
treatment is essentially incorrect, since the orbital order- 
ing is not due to the simple alternation of two arbitrary 
kinds of orbitals, as shown in the following. 





a-orbital 


b-orbital 







3z^ - 


7r/3 


2,y'' - 


z^ — X'' 


27r/3 


i/ - z-^ 


Sx"' — 


TT 


3z^ - 




47r/3 


— X'' 




57r/3 


3x^ - 





TABLE I. Phase and the corresponding eg-electron 
orbitals. 

By the above transformation, H^^^^^ and i?ci-ci are, 
respectively, rewritten as 

^{Sbr(-2(ni>ni + ( 



rrMF 



SjT[2gi(nia - nib) + 9?]}, 



and 



i (iJ) 



(20) 



(21) 



where fii-^ 
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and fii — nia+n-ib- Note that 
i/ci-oi is invariant with respect to the choice of ^i. 
Now let us apply the decoupling procedure as niariib ~ 
(nia)nib + nia(nib) - ("-ia)("-ib), by noting the relations 
("ia) = ((n-i) - gi)/2, (nib) = ((ni) -I- qi)/2, and ni = fii. 
Then, the electron-electron interaction term is given in 
the MFA as 
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H^E,, = {U'/4)Y}2{ni)hi~{ni)^ + 2qi{nai-hM) + q?] 

i 

+ F^[(ni+a)ni - (l/2)(ni+a)(ni)]. (22) 

ia 

By combining H^^^^^ with H^]:^^ and transforming Cia 
and Cib into the original operators as Cia and Cib, the 
mean-field Hamiltonian is finally obtained as 

^UF = - X! ^77''^i7Ci+a7' + -^AF ^ Si • Sj 
ia77' (ij) 

+ -BjT 5Ih2((r.i)T,i + (r,i)T,i) + (r.i)^ + (r,;)^] 

i 

+ J2[iU'/2){n,)+Vj2{n,+^)]n, 

i a 

- 5^[(t/74)(ni) + (F/2)^(ni+a)](ni), (23) 

i a 

where the renormalized JT energy is given by 

i^jT = EjT + U'/A, (24) 

and the renormalized inter-orbital Coulomb interaction 
is expressed as 

= 4£;b,.. (25) 

Physically, the former relation indicates that the JT en- 
ergy is effectively enhanced by U' . Namely, the strong 
on-site Coulombic correlation plays the same role as that 
of the JT phonon, at least at the mean-field level, indicat- 
ing that it is not necessary to include U' explicitly in the 
models, as has been emphasized by the present authors 
in several publications. The latter equation for U' means 
that the one-site inter-orbital Coulomb interaction is ef- 
fectively reduced by the breathing-mode phonon, since 
the optical-mode phonon provides an effective attraction 
between electrons. The expected positive value of U' in- 
dicates that eg electrons dislike double occupancy at the 
site, since the energy loss is proportional to the average 
local electron number in the mean- field argument. Thus, 
to exploit the gain due to the static JT energy and avoid 
the loss due to the on-site repulsion, an Cg electron will 
singly occupy the site. 

If /3 is unrealistically small and i?br becomes much 
larger than [/', U' can be negative, leading to an effec- 
tive attraction between eg electrons, and charge-density- 
wave (CDW) state or s-wave superconducting phase will 
appear. If the dynamical effect is correctly taken into 
account beyond the adiabatic approximation, the com- 
petition among CDW, spin-density- wave, and d-wave su- 
perconducting states will occur .t3 This is an interesting 
point, but in thej-actual manganite /3 is suggested to be 
larger than unity.E3 It has been shown that the breathing- 
mode distortion is suppressed for a reasonable choice pi 
parameters, even if it is included in the calculation.E^I 
Thus, in the following, /? is taken to be infinity for sim- 
plicity, and the breathing mode is dropped from the cal- 
culation. 



B. JT-phononic vs. Coulombic scenario 



In order to check the validity of the present mean-field 
treatment, it is necessary to compare the mean- field re- 
sult with some exact solutions. For this purpose, numeri- 
cal techniques are here applied in small-size clusters. For 
a set of lattice distortions, {gn, g2i, 93i}, the ground-state 
energy is calculated by using the ED method to include 
exactly the effect of electron correlations. By searching 
for the minimum energy, the set of optimal distortions 
{^if* 1 921^*1 'i'sf*} is determined. In this exact treatment, 
the phase to characterize the local orbital set is defined 
by t&n-^iqlT Iql?^). 

In this subsection, a small 4-site ID chain with the 
periodic boundary condition is used in order to grasp 
the most basic aspects of the problem with a minimum 
requirement of CPU time. The small cluster size and 
the low dimensionality will be severe tests for the MFA. 
However, if the validity is verified under such severe con- 
ditions, the mean-field result can be more easily accepted 
in the large cluster limit and in higher dimensions. In or- 
der to focus our attention on the roles of JT phonon and 
on-site correlations, the t2g spin configuration is assumed 
to be ferromagnetic (thus, Jaf will not play an impor- 
tant role) and V is set as zero. The effect of V will be 
discussed in the next subsection in the context of the CO 
state of half-doped manganites. The realistic hopping 
matrix set is used for the ID chain in the x-direction, al- 
though the conclusions are independent of the direction 
of the ID chain. 

In Figs. 2(a)-(f), the results for n=\ and EjT=t are 
shown. Figure 2(a) contains the ground-state energy 
plotted as a function of U' for a fixed value of i?jT . In 
general, the many-body effects cannot be perfectly in- 
cluded within the MFA. However, when U' is introduced 
to the ID chain with Ej^y^O, the results becomes exact 
due to the special properties of one diaiension and the use 
of the realistic hopping amplitudes Namely, the hop- 
ping direction is restricted only along one axis, and the 
orbital set is uniquely determined due to the optimization 
of the JT distortion. In such an orbital polarized situa- 
tion, the electron correlation is included exactly even in 
the mean-field level. Unfortunately, for 2D or 3D FM 
phases, there is no guarantee that the MFA provides the 
exact results, as discussed briefly in the final paragraph, 
but the MFA for the Coulomb interactions is still ex- 
pected to provide the qualitatively correct tendency due 
to the general expectation that the MFA becomes valid in 
higher dimensions. However, in the purely JT-phononic 
model H"^, a remarkable fact appears. Namely, the MFA 
is always exact irrespective of the cluster size, electron 
density, and dimensionality in H^. This is quite nat- 
ural, since the static distortion gives only the potential 
energy for the eg electrons, and this fermionic sector is es- 
sentially a one-body problem, even though the potential 
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should be optimized self-consistently. The lattice distor- 
tion is basically determined by the local electron density, 
and the MFA becomes exact in the static limit for the 
distortion. 




0.0 0.0 

02468 10 012345 

U' (units of t) Ej.^ (units of t) 



FIG. 2. Mean-field and exact-diagonalization results at 
n=l for a 4-site ID chain with the realistic hopping ampli- 
tude along the x-axis. In all figures, the curves denote the 
mean-field results and open symbols indicate the results ob- 
tained by the combination of the exact diagonalization and 
the relaxation method in the phonons. Figures 2(a)-(c) indi- 
cate the (y'-dependence of some quantities for the special case 
E}T=t, while these same quantites are plotted as a function of 
E}T for U'=5t in Figs. 2(d)-(f). In (a) and (d), ground-state 
energies are shown. In (b) and (e), orbital densities, (nia) 
and (nib), are plotted. The solid curve and open square de- 
note (riia) for site 1. The broken curve and open circle denote 
(riia) for site 2. The dash-broken curve and open diamond 
indicate (nib) for site 1. The long-dashed curve and open tri- 
angle indicate (nib) for site 2. In (c) and (f), the phases are 
plotted. The solid line and open square denote for site 1. 
The broken curve and open circle denote for site 2. Note 
that only the results at site 1 and 2 are depicted, since the 
results at site 3 and 4 are the same as those at site 1 and 2, 
respectively. 

For the reasons discussed above the MFA works quite 
well in the JT scenario, but this fact does not mean that 
the obtained state is trivial, since the orbital degree of 
freedom is active and several non-trivial 00 phases oc- 
cur. To visualize this result, the orbital densities, (riia) 
and (nib), are shown in Fig. 2(b), and the phase ^i is plot- 
ted in Fig. 2(c). It is observed that (nia) becomes very 
small and (nib) is almost unity for all the sites. For even 
and odd sites, ^i is given by 27r/3 and Stt/S, respectively. 
Namely, the occupied b-orbital is d^^2_j,2 at the even sites 
and dy2_^2 at the odd sites, respectively. This is simply 
the orbital-staggered state in the spin FM phase, which 



has been observed in the MC analysis, til The mechanism 
of its occurrence is quite simple: Imagine the two adja- 
cent sites in which one Cg electron is present at each site. 
If the limit Ej^'^t is considered, the energy gain can be 
evaluated in second order with respect to and this en- 
ergy is maximized when the occupied orbital in a certain 
site is exactly the same as the unoccupied orbital in the 
neighboring site. Namely, the difference in ^i between 
two adjacent sites should be tt at n—1. Note that the 
value of is determined by the hopping direction. 

Next let us examine the Coulombic scenario, still us- 
ing a small 4-site cluster to compare mean-field results 
against exact ones. In Figs. 2(d)-(f), the ground-state 
energy, orbital densities, and the phase are plotted as a 
function of Ejt for a fixed value of U'=5t. Again it can 
be observed that the MFA works quite well, except for 
the case of £'jt=0. In this pure Coulombic limit, the 
JT distortion docs not occur, and thus, the results are 
obtained only by the ED for H^. The ground state en- 
ergy obtained in the MFA converges to the exact result 
in the limit of E^t^ 0. However, no symbols are shown 
at £'jT=0 in Figs. 2(e) and (f), since the orbital densities 
and the phase could not be fixed at £'jt=0. This is due 
to the fact that the energy is invariant for any choice of 
local orbital, indicating that (n-ia), (nib), and ^i cannot be 
determined at each site for £'jt=0. In other words, the 
00 state is not uniquely fixed within the purely Coulom- 
bic model in the sense that the a special orbital is not 
specified at each site. However, the orbital staggered ten- 
dency is detected if the orbital correlation as a function 
of distance is studied, as will be discussed in Sec. VI. 

Now the roles of the JT-phonon and on-site correlation 
are discussed. At finite electron-phonon coupling, the op- 
timized orbital is determined, and the MFA provides an 
essentially exact result for the shape of the orbital. Note 
that carrying out unbiased MC simulations and using the 
relaxation technique are still important tasks, since the 
present MFA works quite well only for a fixed t2g spin 
background. Thus, it is an unavoidable step to check 
whether the assumed t2g spin pattern is really stable or 
not with the use of unbiased techniques by optimizing 
the lattice distortions as well as the t2g spin directions 
simultaneously. In this sense, the analytic mean-field ap- 
proach becomes very powerful when it is combined with 
numerical unbiased techniques. 

In the purely Coulombic model, however, no special 
orbital is determined at each site due to the local SU(2) 
symmetry. Only when the JT distortion is included, the 
optimal orbital set at each site is determined from the 
competition between the kinetic energy gain and the po- 
tential loss due to the lattice distortion. As a conse- 
quence, it appears that the natural starting model to 
understand the properties of manganites should be the 
JT-phononic model, and it is enough to include the ef- 
fect of U' on Hj^. Although the validity of this statement 
has been shown in a limited situation in this paper, it is 
believed to be correct in other cases. (As for the metal- 
insulator transition, a comment will be given in the final 
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paragraph of this subsection.) On the other hand, if the 
starting model is chosen as the purely Coulombic model, 
the degeneracy in the orbital space is lifted when the JT 
distortion is introduced, indicating that the ground-state 
property abruptly changes due to the inclusion of the JT 
distortion. However, the on-site correlation is still im- 
portant to achieve effectively the strong-coupling region, 
since U' is renormalized into the JT energy. 




2 4 6 
U' (units of t) 



2 3 4 
(units of t) 



FIG. 3. Mean-field and exact-diagonalization results at 
n—0.5 for a 4-site ID chain with the realistic hopping am- 
pfitude along the a;-axis. The meaning of curves and symbols 
are the same as those in Fig. 2. 

To check whether the statements in the previous para- 
graph are correct also in the doped case, the same anal- 
ysis is carried out for n=0.5, as shown in Figs. 3(a)-(f). 
Note again that the MFA results are always exact and 
the physical quantities are continuous as a function of U' 
for non-zero values of Ejt, but the local orbital set can- 
not be uniquely determined for i?jT=0. In addition, the 
results are independent of U' in Figs. 3(a)-(c). This is un- 
derstood as follows: If a unitary transformation is intro- 
duced to make the hopping matrix diagonal, the model 
Hq' in the FM-phase is reduced to 



Hfk - -^(4^/3)a^ai+a- 



C/'^a^ai6|5i, 



(26) 



where Ui and bi are given by ai=—{^/3/2)cia+{l/2)cib 
and 6i=(l/2)cia-l-(\/3/2)ci6, respectively. In this reduced 
model, only the "a" -electrons can hop and they interact 
with the localized, "b" -electrons, which is just the Falicov- 
Kimball model.E^ Of course, this model is obtained triv- 
ially if the ID chain along the z-direction is considered. 
Thus, the ground state at n=0.5 is obtained by filling the 
lower "a"-band up to 7i=0.5 and U' does not work at all. 



When the electron-phonon coupling is switched on, the 
system becomes insulating, but the effect of U' is still 
inactive. Although the perfect agreement of the mean- 
field results with the exact numbers is due to the special 
properties of the ID chain, the JT-phononic model with 
the realistic hopping appears to contain the important 
physics of the problem and the electron correlation is 
not as crucially important as naively expected. 

As for the orbital densities at n=0.5, (nia) is always 
negligible small, but (nib) for the even sites is almost 
unity, while (nbi) for the odd sites is small, clearly indi- 
cating the CO tendency. In this case, is always given 
by 27r/3 irrespective of the site position, denoting the 
occurrence of the ferro ordering of the d^^^^r^ orbitals. 
This is easily understood by the DE mechanism in the 
orbital degree of freedom. Namely, the orbital arranges 
uniformly to improve the kinetic energy of the eg elec- 
trons, just as the spin does in the FM phase of the doped 
one-orbital model. 

Finally, a brief comment on the MFA in higher di- 
mensions is provided. Fop-i7'=0, the ID system is insu- 
lating if i?jT is non-zero,C3 but in higher dimensions, a 
metal-insulator transition occurs at some finite value of 
EjT. If EjT is larger than this critical value, the ground 
state wave function will be smoothly changed as U' in- 
creases. However, when Ejt is so small that the sys- 
tem is metallic at U'—O, the inclusion of U' will bring 
the metal-insulator transition at some value of U'. The 
MFA cannot predict correctly this metal-insulator tran- 
sition point, but in this paper, the CO state is mainly 
discussed. Namely, U' (Ejt) is included into the insulat- 
ing phase originating from Ej^iU'). Although the prop- 
erties of the metal-insulator transition in the manganese 
oxide are quite interesting, this type of analysis will be 
left for the future. 



IV. CHARGE-ORBITAL ORDERING IN 
HALF-DOPED MANGANITES 

In the previous section, the CO-state with a uniform 
orbital arrangement has been suggested as the ground- 
state of the ID chain at n=0.5. However, in the real 
materials a more complicated situation occurs. In the 
narrow-band manganite such as Lao.sCao.sMnOs and 
Ndo.sSro.sMnOs, the CE-type AFM phase is stabilized 
(See Fig. 4(a)). In this structure, the t2g spins form 
a ferromagnetic array along a zigzag ID path, the CO 
state appears with the checkerboard pattern, and the 
(3x^ — r^)/(3?/^ — r^)-type orbital ordering is concomitant 
to this CO state. The origin of this complex spin-charge- 
orbital structure has been recently clacified on the basis of 
the topology of the zigzag structure.E3 Along the z-axis, 
the Ci?-pattern stacks with the same CO/00 structure, 
but the coupling of i2g spins along the z-axis is anti- 
ferromagnetic. If only the bi-partite charge structure in 
the x-y plane is considered, naively the nearest-neighbor 
Coulomb interaction V seems to be the key issue in the 
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formation of the CO/00 phase.Ea However, if V is strong 
enough to bring the CO structure in the x-y plane, the 
3D Wigner-crystal structure should appear in the cubic 
lattice, but this is not observed in the real system. Thus, 
V is not the only ingredient needed for the stabilization 
of the CO structure of half-doped manganites, as already 
discussed in the Introduction. What other terms in the 
Hamiltonian may create a charge-stacked phase? 



(a) 



■• * ■ 
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FIG. 4. (a) Schematic representation of the spin, charge, 
and orbital ordering for the C-E-type AFM structure. Solid 
and open symbols indicate up and down t2g spins, respec- 
tively. The lobes indicates {Sx^ — r^)- or (3j/'^ — r^)-orbitals 
at the Mn^+ sites, while the circles denote Mn*''" sites, (b) 
Energy per site as a function of Jaf for A=1.6 and U'=0. 
The curves denote the mean-field results and the solid sym- 
bols indicate the energy obtained by the relaxation method. 
Thick solid, thick broken, thin broken, thick dashed, thin 
dashed, thin broken, and thin solid lines denotes FM, A-type, 
shifted C-B-type, C-E-type, C-type, and G-type states, respec- 
tively, (c) Phase diagram in the {Jaf,V) plane. Note that 
the charge-stacked structure along the 2-axis can be observed 
only in the C-B-type AFM phase. 

In order to clarify this point, in this subsection 4x4x4 
lattices with the periodic boundary condition are stud- 
ied as a simple representation of half-doped manganites 
on the basis of H°°. In Fig. 4(b), energies for several 
magnetic phases are plotted as a function of Jaf for 



A = T6B At this stage in the calculations, V is set 
to zero and its effect will be discussed later. The results 
of Fig. 4(b) are obtained from the JT-phononic model, 
but as mentioned in the previous section, it can be inter- 
preted that the effect of U' is included effectively in the 
electron-phonon coupling. Although the effect of U' also 
appears through the non-trivial term (C/'/2) ^;(ni)ni, 
this term essentially indicates the prohibition of double 
occupancy, and does not lead to qualitative changes in 
the results obtained from the JT-phononic model, as long 
as /3 is larger than unity and U' is positive. Thus, the re- 
sults below can be interpreted as arising from a purely JT 
calculation or a mixture of JT and Coulomb interactions 
within a mean-field technique. The curves shown are the 
results in the MFA obtained for several fixed t2g spin 
patterns, while the solid circles are obtained by the opti- 
mization of both the local JT distortion and the t2g spin 
angles simultaneously. The agreement between the an- 
alytic and numerical results is excellent, indicating that 
the MFA works quite well in the JT-phononic model. For 
JAF^O.li, a metallic 3D FM-phase is stabilized, since 
this spin pattern optimizes the kinetic energy. In a very 
narrow region around JAF~0.lt, the ^-type AFM phase 
occurs. This phase is metallic in the present intermediate 
coupling, but the uniform (x^ — y^)-type orbital ordering 
appears. In the wider region O.lt^ JAF^0.25t, the C-E- 
type AFM structure is the ground state. For unrealistic 
large values of Jaf, the G-type AFM phase is stabilized 
to gain the magnetic energy. 

Now let us focus our attention on the stabilization of 
the C-E-type phase. Due to the competition between 
the kinetic energy of the Cg electrons and the magnetic 
energy of the t2g spins, the, one-dimensional stripe-like 
AFM configuration occurs ,Ej in which arrays of t2g spins 
order ferromagnetically along some particular ID paths. 
Precisely the shape of the ID FM path, the zigzag path, 
is a key to understand the CE-type phase. To see this, 
let us compare the CE-state with the C-type AFM state, 
which is characterized by the straight-line FM path. Al- 
though the energy of the C-type structure has the same 
slope as a function of Jaf as the CE-type, it is not the 
ground-state, as shown in Fig. 4(b). As for the orbital ar- 
rangement, the (3a;^ — r^)-type 00 occurs in the C-type, 



while the (3x 



■)-type orbital arrangement 



appears in the C-E-type state. At a first glance, these two 
00-states seem to be quite different, but if the concept 
of the orbital DE mechanism is employed, it is reason- 
able that there is no essential distinction between them. 
Namely, in both cases, the Cg electron orbital is always 
polarized along the hopping direction of the ID paths. 

An essential point for the stabilization of the charge- 
stacked CE-type AFM phase is the difference of the 
topology of the ID paths between the straight and zigzag 
shapes. Mathematically, the topology of the IDjjath can 
be characterized by "the winding number" itllj defined 
from the Berry phase connectionpOf the Cg-electron wave 
function along the hopping pathMd As shown in Ref. |2^, 
w is decomposed into two terms as w^Wg+Wt- The for- 
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mer, Wg, is called the geometric termJ13 which is (1) 
corresponding to the ferro- (antiferro-) arrangement in 
the orbital sector along the ID path. As discussed above, 
due to the orbital DE mechanism, Wg is zero in the doped 
case irrespective of the shape of the hopping path. The 
difference between the straight and zigzag paths appears 
in Wt, expressed as wt=N^/2, where is the Jiumber 
of vertices appearing in the unit of the ID chain.tJ Since 
Wt is determined only by the shape of the ID path, it is 
called the topological term. Of course, A'v is zero for the 
straight path, while is equal to 2 for the zigzag path. 
Thus, it is obtained that 'w=Q and 1 are the values for 
the C- and Ci?-type AFM phases, respectively. 

To understand why the zigzag chain with w—\ has 
the lower energy, it is useful to consider the situation 
U'=E}'T={), in which the Ci?-type AFM spin structure 
characterized by the zigzag ID cha in is stabilized in the 
picture of a band insulator. E2Ie§cj Due to the periodic 
change along the zigzag path in the hopping amplitude 
as { • • •, i^^,, t^^,, t^^,, t^y, ■■■ }, the system becomes 
the band insulator. Although the difference in the hop- 
ping amplitudes in the x- and y-directions is only a phase 
factor in t^^^ and t^j^, a bandgap of the order of t de- 
velops. On the other hand, note that the C-type AFM 
states characterized by the straight path is metallic, since 
there is no change in the hopping amplitude from bond to 
bond. In Ref. PQ, it has been clearly shown that the CE- 
type structure has the largest bandgap among all the pos- 
sible types of zigzag hopping paths at n=0.5. Moreover, 
in Ref. if the JT energy is switched on in this phase, it 
has been shown that the Ci^-type phase continues to be 
the ground state and charge-ordering appears with the 
{3x^ — r^)/(3?/^ — r^)-type orbital-ordering. Thus, based 
on the band-insulator picture and the spirit of the adia- 
batic continuation, the Ci?-type AFM structure charac- 
terized by the ID zigzag FM chain is the ground-state. 

Now the physical meaning of the energy difference be- 
tween CE- and C-type states is discussed using the the 
concept of the winding number. For this purpose, an 
analogy with a typical spin_4j)roblem is quite useful, as 
suggested by Takada et alE3 In the spin problem, by 
classifying the states with the total spin 5" which is a 
conserved quantity, the exchange energy is defined by 
the difference between the energies of the singlet {S=Q) 
and triplet (>5'=1) states. In the present problem, w is 
the topological entity and the conserved quantity. Thus, 
the energy difference between the states with w — 1 and 
w = 0, Jw, is expected to play a similar role as the ex- 
change energy in the spin problem. As for the magnitude 
of Jw, it can be of the-|Order of O.lt from the analysis of 
the two-site problem,E3 although it depends on the value 

of E]T- 

Another phase competing with the £2£Ltype state is 
the shifted CE-type {sCE) structurcH-El This is also 
obtained by the stacking of Fig. 4(a) along the 2;-axis, 
but one-lattice spacing shifted in the x- (or ?/-)direction. 
Due to this shift, the number of FM and AFM bonds 
becomes equal, and the magnetic energy is exactly can- 



celed. In fact, the energy for the sCE structure is inde- 
pendent of Jaf- The Ci?- type phase is stabilized against 
the sC-B-structure by the magnetic energy, as observed in 
Fig. 4(b), showing the key role played by Jaf in models 
for manganites and likely in the real compounds as well. 

Let us consider now the effect of V on the CO states. 
For this purpose, the mean-field calculation is carried out 
based oii|-ff^p for a reascmable parametei seta such as 
t=0.5eV,N £;jT=0.25eV,El and {7'=5eV.El The phase di- 
agram in the {Jaf, V) plane is shown in Fig. 4(c). From 
this figure it is clear that the CS structure occurs for 
the CiJ-type AFM phase in a broad region of parame- 
ter space, as deduced from the Fourier transform of the 
charge correlation which has an observed peak at (tt, tt, 0) 
in the C-E-type state, while a peak-appears at (7r,7r,7r) 
in the sCE and C-type AFM statesO A remarkable fact 
of Fig. 4(c) is that the CS structure is robust against 
the inclusion of V. The boundary between the sCE- 
and Ci?-phases is qualitatively understood by the bal- 
ance of tjbe magnetic energy gain and the charge repul- 
sion IossBj On the other hand, the phase boundary be- 
tween the CE- and C-type AFM states is independent 
of Jaf, since those two states have the same magnetic 
energy. As mentioned above, in this case, the energy J„ 
due to the difference in the topology of ID path stabi- 
lizes the CS structure in spite of the charge repulsion 
loss. In fact, the phase boundary exists around Vf^Q.3t, 
which is the same order as Jw. Although it is difficult to 
know the exact value of V in the actual material, if the 
simple screened Coulomb interaction is estiHiated using 
the large dielectric constant of manganites, eZI V/t is esti- 
mated to be 0.1~0.2,B i.e., inside the CS region m our 
phase diagram. 

Finally, a comment on the stabilization of the Ci?- type 
structure in the purely Coulombic model is provided. As 
discussed above, in the case of /7'=i?jT=0, the Ci5-type 
AFM spin structure is found to be the ground-state in the 
band-insulator picture. If U' is smoothly switched-on in 
this phase, still keeping Ejt—0, the Ci?-type phase still 
continues to be the ground state and charge-ordering ap- 
pears even without the help of V, since the local charge 
in the straight segmeni_pf the zigzag chain is larger than 
that at the corner site.EJ The ground state energy for the 
zigzag chain is lower than that for the straight chain, and 
its difference is again of the order of O.lt. Since this en- 
ergy difference can compensate the energy loss due to V 
in the CS structure, it would be possible to understand 
the CS structure even in the purely Coulombic model on 
the same topological argument as carried out for the JT- 
phononic model. Thus, it is possible to fill the 3D cubic 
lattice by the zigzag ID chains stacked in the b- and c-axis 
directions, with the same charge ordering but antiparal- 
lel t2g-spin directions across those ID chains. Note, how- 
ever, that the energy is invariant for the choice of local Cg- 
electron orbital and the (Sa;^ — r^)/(3y^ — r^)-type orbital 
arrangement cannot be specified in the purely Coulom- 
bic model. On the other hand, in the pure JT-phononic 
model, the Ci?-type AFM spin structure, the charge- 
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stacked CO state, and the (Sx^ — r^)/(3i/^ 



-type 00- 



phase have been fully understood, as explained before in 
this section. Thus, based on these results, these authors 
believe that the purely JT-phononic model is more effec- 
tive than the purely Coulombic model for the theoretical 
investigation of manganites, although certainly both lead 
to very similar physics. 



V. PHASE-SEPARATION TENDENCY 



In previous investigations using the unbiased MC sim- 
ulations, the PS teadeiiey has been clearly established in 
manganite models .QcIO In the simulatioas- PS appeared 
both in the one-orbital FM Kando modeflij and the two- 
orbital JT-phononic models.llil PS occurs due to the bal- 
ance between the kinetic energy of the Cg electrons and 
the potential energy due to the background, either the 
t2g spins or the JT distortion. Due to this difference in 
the origin of the background potential, two-types of PS 
exist, spin driven and orbital driven. In the following, 
the ID case is used for simplicity. In higher dimensions, 
the situation will be more complicated, but the essential 
physics is expected to be captured in the ID case. 

One type of PS appears in the two-orbital model driven 
by the spin sector between the spin FM and AFM phases, 
mainly in the region between xwO.S and x—1.0. To im- 
prove the kinetic energy, in this doped situation the t2g 
spins and eg-electron orbital tend to array in a ferro- 
manner, but in the heavily doped region close to a;=1.0, 
the t2g spins order antiferromagnetically to gain the 
magnetic energy which dominates over the kinetic en- 
ergy. Thus, in this case, the PS between FM/OF- and 
AFM/OF-states appears, as clearly indicated in the pre- 
vious MC simulations. Note that the acronym in front of 
the slash indicates the spin state, while the acronym after 
the slash denotes the orbital arrangement, e.g., FM/OF 
indicates the spin FM and orbital ferro state. Note that 
this PS is possible in the one-orbital model as well, since 
the orbital degrees of freedom is not active for this case. 

Another variety of PS is related to the orbital degrees 
of freedom, which is tightly coupled to the JT distor- 
tion, which mainly exists between a;=0 and a;«0.5 in the 
two-orbitals model. In the undoped situation and for rea- 
sonable values of Jaf, the FM/OAF-state is the ground 
state, as shown in the previous section. Also in the up^ 
biased MC simulation, this phase has been obtained.tiJ 
When holes are doped, the spin structure is still FM, but 
the orbital arrangement becomes uniform to gain the ki- 
netic energy by polarizing the orbitals along the hopping 
direction (orbital DE mechanism). Thus, in this case, 
the PS between FM/OAF- and FM/OF-states appears, 
as suggested in the previous works for the two-orbital 
model.llil 

In the MC simulation, the grand canonical ensemble 
has been used and the chemical potential /i was tuned to 
obtain the target electron density. In such a calculation, 



the presence of PS was clearly observed by monitoring 
the density vs. fi. There is a range of densities that can- 
not be stabilized, no matter how carefully /i is fine-tuned. 
In other words, the plot (n) as a function of /i presents 
a discontinuity at some particular value of fi. This effect 
occurs for the one- and two-orbital models, the latter 
with JT-phonons, and in all dimensions of interest. Fur- 
ther evidence of mixed-phase tendencies can be obtained 
by monitoring the density as a function of MC-time. In a 
stable regime, the density does not change much between 
MC configurations, but in the PS region the fluctuations 
are very strong. 
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FIG. 5. (a) Energy difference AE as a function of x for a 
16-site ID chain with U'—Q and i5jT=1.2. The definition of 
AE is shown in the text, (b) Orbital density (nib) for x—0 
(upper chain) and 0.5 (lower chain). Note that the b-orbital in 
the x-y plane is depicted and its size is proportional to (nib), 
(c) Energy difference AE as a function of x for a 16-site ID 
chain with U' — W and i5jT=0. 

In the present MFA, the electron number n has been 
fixed. For a given n, the optimized structure both for the 
t2g spins and the lattice distortion are determined by the 
analytic MFA and the numerical technique. To under- 
stand the PS tendency in this formalism, it is necessary 
to check the stability of the obtained phase by calculating 
the ground-state energy Eq as a function of n. Namely, 
if a negative curvature is obtained, i.e., d^EQ/dn^<0, 
such a phase is unstable even if it is the ground state at 
a fixed n. If the PS tendency is included in the present 
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model, the FM/OF states between a;=0 and x=0.5 should 
be unstable and the mixed phase should have lower en- 
ergy than the states obtained by the MFA, in order to 
reproduce the MC results. In order to verify this, the 
energy difference AE{x) is plotted as a function of x for 
Ej'T — 1.2 in a 16-site spin-FM ID-chain with the real- 
istic hopping along the x-axis. Here AE{x)=EQ{x)—ex 
withe=2[£'o(x=0)-£'o(a;=0.5)]. If AE{x) is positive, the 
homogeneous state is unstable and, instead, the mixed- 
phase appears. Note that the state at a;=0.5 is stable. 
The result is shown in Fig. 5(a), in which positive AE{x) 
is indeed observed between x=0 and x=0.5. In Fig. 5(b), 
the orbital arrangements at a;=0 and x=0.5 are shown 
by depicting the shape of the occupied b-orbital in the 
x-y plane. It should be noted that the size of the or- 
bital is proportional to (hbi). This result agrees very well 
with the previous MC calculation, and the PS tendency 
is believed to be definitely established in the JT-phononic 
model. 

Now let us analyze whether it is possible to detect the 
PS tendency in the purely Coulombic model or not. For 
this purpose, A£'(a;) is evaluated using a 16-site spin-FM 
ID-chain with the realistic hopping along the x-axis for 
^'JT = and U' = lOt by using the mean-field Hamilto- 
nian Eq. (p3|). As shown in Fig. 5(c), again the positive 
AE for X between and 0.5 can be observed, indicat- 
ing clearly the PS tendency. If the present mean-field 
Hamiltonian is accepted, this result is quite natural, since 
U' is included effectively in the coupling between the Cg 
electrons and JT distortion. However, one may consider 
that this is just an artifact due to the MFA. Thus, in the 
following section, it is explicitly shown that this PS ten- 
dency in the purely Coulombic model is not due to par- 
ticular properties of the MFA by performing the DMRG 
calculation in the ID Hubbard-like model. 

In short, it is quite interesting to observe that the MFA 
can properly reproduce the PS tendencies found using 
other more sophisticated techniques. Although the cal- 
culation is fairly simple and handy the obtained result is 
physically meaningful and reliable. More investigations 
on the PS tendency by using the MFA in higher dimen- 
sions would be certainly interesting since the MC sim- 
ulations become increasingly difficult as such dimension 
grows, but this point will be discussed in future publica- 
tions. 



VI. DMRG RESULT FOR He 

The purpose of this section is to continue the analy- 
sis of the reduced Hamiltonians of Section II, this time 
focusing on the purely Coulombic model He- The multi- 
orbital Hubbard model has-been addressed before using a 
variety of approximations but here care must be taken 
to select the appropriate technique accurate enough to 
search for physics similar to the results obtained with 
JT-phonons. The use of unbiased methods is particularly 



important for subtle issues such as PS. For this purpose, 
the model Hq described in Siection II has been studied 
here with the DMRG method,E3 supplemented by the ED 
technique, keeping the truncation errors around 10~® by 
using typically 120 states on intermediate size chains with 
open boundary conditions. The restriction to work in the 
ID system is not severe in view of the results of previous 
sections and those of Refs. that showed a strong sim- 
ilarity between one, two and three dimensions, at least 
regarding the rough features of the ground state. Actu- 
ally, it is important to remark that the phase diagrams 
described in previous work by our group are mostly based 
on evidence coming from the short-distance correlations 
calculated in our study which are not expected to depend 
strongly on the dimensionality. In ID systems, the ex- 
istence of genuine long-range order vs. slow power-law 
decay of correlations is a subtle issue beyond the goals of 
the present analysis. 

In this context, the static observables studied with the 
DMRG method are the spin structure factor 

S{k)^^Y.(^,-s^)e^(^'^^\ (27) 

the charge structure factor 

^(fc) = ^E(/'^-^'")e'^'""^'' (28) 
and the orbital structure factor 

where k is momentum, j and m denote site positions, 
L is the length of the ID chain, Tj={pja—pjh)/2, and 
Pj-f=Ylia Pjii^- important to remark that the anal- 

ysis described below has focussed on the analogies with 
the results found using JT-phonons, especially regard- 
ing PS and the existence and properties of the spin-AF 
orbital-staggered phase at density n=l. In the several 
other interesting phases reported in this section, our ef- 
fort has been limited to the description of their main 
features regarding their spin, charge and orbital char- 
acteristics, postponing for a future publication a more 
detailed analysis of its origin and possible relevance to 
experiments. In order to focus on the similarity in the 
effects of the on-site correlation and the JT phonons, V 
is set to be zero for the time being. The influence of V 
will be discussed later. Note that the energy unit is t in 
this section, as in the previous ones. 

A. Unit Hopping Matrix 

Thus far the realistic hopping matrix with non-zero off- 
diagonal hopping amplitude has been used. This hopping 
matrix is quite important to understand the properties 
of manganites, but the analysis is sometimes complicated 
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due to the non-zero ofF-diagonal element. In order to an- 
alyze the two-orbital model more easily, it is useful to 
introduce first a simple unit hopping matrix for any di- 
rection, given by 

t^y=tS^Y, (30) 

where Sij is the Kronecker's delta. Besides its simplicity, 
this hopping amplitude has been used before in the search 
for the FM state in the nmlti-orbital Hubbard model with 
a variety of techniques,Ej and thus, it is meaningful to 
investigate its properties. In addition, it will be the- 
oretically interesting to compare results using different 
hopping sets, to study the dependence of the ground- 
state properties with those amplitudes. It will actually 
be concluded that the use of realistic hoppings is very 
important in this context. 




1. Case of n—1 

Consider first the special case of density n—1. Us- 
ing the DMRG and ED methods a large set of cou- 
plings {U', J) have been investigated, but here only the 
most representative results are presented as examples. 
In Fig. 6(a), the spin structure factor S{k) is shown at 
fixed U'—IQ, for particular values of J, and three regimes 
are clearly identified. At small J compared with U' , the 
spin sector has incommensurate characteristics, with a 
maximum at momentum k—TT/2. As J grows, an abrupt 
transition to the FM state is observed, with S{k) now 
peaked at zero momentum. This last regime occurs in a 
robust window of J. At J~15, the spin sector becomes 
incommensurate again with a peak in S{k) at /c=7r/2. 
The charge structure factor N{k) is shown in Fig. 6(b) 
for the same set of parameters. It is observed that at J— 2 
and 5 the charge correlations are not enhanced since only 
a broad peak with low-intensity appears at k=Tr. How- 
ever, for J— 12 and 20 indications of charge-ordering ten- 
dencies are found, since now the k—i: result is clearly 
enhanced, indicating a structure with an alternation be- 
tween the even- and odd-sites of the chain. Note that 
J=5 and 12 have different characteristics when analyzed 
using iV(fc), while they are both ferromagnetic according 
to S{k). 

In Fig. 7(a), the orbital structure factor T^{k) is pre- 
sented for the same set of couplings as used in Figs. 6(a) 
and (b). In this case, enhanced orbital correlations exist 
in the ground state at J— 2 due to the robust values that 
T^{k) has, particularly at k—n. This is indicative of a 
staggered arrangement of orbitals, similar to the results 
discussed before in Sec. III. A qualitatively similar but 
even more prominent effect occurs at J=5. At J=12 and 
20, T^{tt) is relatively small and there is no noticeable 
indication of enhanced orbital correlations. The transi- 
tion from the low-J regime, exemplified by J=2, to the 
intermediate one ( J=5) is abrupt, as already observed in 
the study of the spin structure factor. 



FIG. 6. DMRG results for (a) spin structure factor S{k) 
and (b) charge structure factor N{k) as a function of k for 
He with the unit hopping matrix in the 20-site ID chain for 
U'=10. The dashed, solid, dotted, and dotted-dashed lines 
correspond to J=2, 5, 12 and 20, respectively. 




FIG. 7. (a) Orbital structure factor r^(fc) vs. k us- 
ing the same parameters as in Fig. 6. (b) Phase di- 
agram of He for the unit-matrix hopping amplitude at 
n=l. The four phases are IC/OAF, FM/OAF, FM/CO, 
and IC/CO. Here "IC" means spin-incommensurate, "FM" 
means spin-ferromagnetic, "OAF" denotes orbital-staggered, 
and "CO" indicates charge- ordred. 

Repeating a similar analysis for several other values 
of U' and J allowed us to sketch a phase diagram for 
the two-orbital Hubbard model at density n=l, which is 
shown in Fig. 7(b). Since the complexity of the Hamil- 
tonian does not allow us to perform a careful finite-size 
study to determine the dominant correlations in the bulk 
ground-state, the phase diagram of Fig. 7(b) should be 
considered only qualitative. Nevertheless, since for clus- 
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ter sizes smaller than L=20 the results were found to be 
similar, no large size-effects are anticipated. The four 
identified regions are marked. Note that the emphasis 
has been given to the determination of the boundary of 
the spin-ferromagnetic orbital-ordered phase due to its 
similarities with the analogous phase reported in the. MC 
studies of the two-orbital model with JT phonons,llil and 
in the MFA of the previous section. The boundaries of 
the other regimes are less accurately determined, partic- 
ularly at small values of J and U' . 

Let us consider the origin of the many phases observed 
in Fig. 7. First address the CO regime which appears 
above the line J=U' . The intuitive explanation for this 
behavior is simple. The inter-orbital exchange interac- 
tion J actually provides an attractive interaction to the 
electrons, which try to doubly populate half the sites of 
the chains when this interaction dominates. On the other 
hand, the inter-orbital repulsion U' prevents double oc- 
cupancy of two orbitals at the same site. Then, a compe- 
tition between the two occurs naturally. If J>U', which 
is an unphysical limit, then charge-ordering occurs with 
roughly half the sites having two particles and the other 
half none. To allow for some nonzero electronic kinetic 
energy, the doubly occupied sites do not cluster together 
but are spread on the chain. The spin and charge pat- 
tern compatible with the results of Figs. 6 and 7 in the 
CO regime has a unit cell of size four lattice spacings, 
and an arrangement (2u,0,2d,0) where "2u" denotes two 
spins up, "0" is an empty site, and "2d" are two spins 
down. Spin staggered patterns appear frequently in the 
CO states. 

The region U'>J is physically more interesting, and 
connected with the results observed for the model with 
JT-phonons at the same density. The spin FM charac- 
teristic of the FM/OAF phase is beheved to be caused by 
the optimization of the kinetic energy by spin alignment. 
The influence of the attractive coupling J is also very 
important for the stabilization of this phase, since it also 
favors spin alignment at every site. The orbital-staggered 
pattern allows for some mobility of the electrons from site 
to site, while an orbital-uniform arrangement would not 
allow for that movement at density n=l, if all spins were 
aligned and when the unit-matrix hopping is used. The 
charge is spread uniformly, i.e., no charge-ordering exists 
in this phase. However, below but close to the line J=U', 
charge correlations are enhanced due to the proximity to 
the CO-regime. In this subregion, the orbital correla- 
tions are suppressed compared with the results observed 
at the lower J end of the FM/OAF phase, where they 
are maximized. Overall, it is clear that the FM/OAF 
state has characteristics very similar to those observed 
when JT-phonons are used to mediate the interaction be- 
tween electrons (See Sec. III). Then, this phase appears 
prominently both in studies with phonons and Coulombic 
interactions, and likely it will be stable when a mixture 
of the two terms is used. 

As J is decreased for fixed U', the ferromagnetic ten- 
dencies are naturally also reduced. The computational 



study shows that a regime with sharp spin incommen- 
surate characteristics dominates in this region. The or- 
bital order remains staggered and the charge is uniformly 
spread. This state is not directly related with the goals 
of the paper, and thus, the discussion of its origin and 
characteristics is postponed for future work. 



2. Case ofn^l 

One of the main goals of the study in this section is the 
investigation of whether PS tendencies appear in a multi- 
orbital model having only Coulomb interactions, starting 
at n=l in the FM/OAF regime previously observed with 
JT-phonons. For this purpose, here the couplings were 
fixed to U'=30 and J=22, i.e., inside the FM/OAF phase 
of Fig. 7(b), and the density was varied using a cluster 
of 20 sites. As the number of electrons N was changed 
between 6 and 28 (only using an even number) , it was ob- 
served that the ferromagnetic characteristics persist and 
S{k) continued to be sharply peaked at zero momentum. 




FIG. 8. (a) DMRG result for the charge structure factor 
N{k) vs. k at n—1 for He with the unit matrix hopping ampli- 
tude, and a 20-site ID chain. U' and J are fixed as 30 and 22, 
respectively, corresponding to couplings inside the FM/OAF 
phase in Fig. 7(b) . The number of electrons A'' is shown in 
the figure. Note that the fc=7r response is enhanced at density 
n=0.5, at least relatively to the other densities, (b) Orbital 
structure factor T^(fc) vs. k using the same convention and 
couplings as in (a). Note the appearance of incommensurate 
characteristics in this channel as the density changes away 
from 71=1. 

Regarding charge-ordering, an enhancement in this 
channel was found at density n=0.5 as can be observed in 
Fig. 8(a) where results for A^=8, 10, 16, and 20 are pre- 
sented. This enhancement shows that at the couplings 
studied here tendencies toward a CO-state are devel- 
oped at 71=0.5, in agreenient with recent MC studies for 
cooperative JT-phonons,E3 the analysis of the previous 
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sections with non-coopjerative phonons, and with exper- 
iments for manganitesE The results in the orbital sector 
are shown in Fig. 8(b). As the density is reduced, or in- 
creased, starting at n—1 the orbital correlations develop 
incommensurate characteristics, which is a curious effect 
not observed before to the best of our knowledge. When 
the density reaches n — 0.5, T^{k) peaks at effect 
likely correlated with the precursors of charge-ordering 
found in N{k). 
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FIG. 9. (a) Density n vs. \i constructed from the 
ground-state energies corresponding to several number of elec- 
trons on a chain with 20 sites, and the couplings used in 
Fig. 7. Note that all densities are accessible for some val- 
ues of [s,. (b) Phase diagram of He for the case of the unit 
hopping matrix at n=0.5. The two phases identified are 
IC/OIC and FM/OIC, where "OIC" means charge-disordered 
orbital-incommensurate. The charge-ordered regime above 
J=U' has not been studied in as much detail as in Fig. 7(b), 
and it is simply denoted by CO. 



Note that with the DMRG and ED techniques, which 
are setup in the canonical ensemble where the number 
of particles can be fixed arbitrarily, the stability of the 
various states with A^-electrons cannot be addressed di- 
rectly. For this purpose it is necessary to compare the 
energies of the various ground states and construct a plot 
of the density n vs. /i following steps already described 
in detail in previous publicationsja and in the previous 
section. Figure 9(a) shows that all the densities studied 
here are actually stable in the sense that a finite window 
of (X exists for all of them where the state under study 
minimizes the energy. This has to be contrasted with 
the results found in the MC and mean-field calculations 
with JT-phonons where states in a finite window of n 
were found to be unstable,tll i.e. there was no value of 
jjL that render them the global ground-state of the sys- 
tem. Then, it is concluded that the two-orbital Hub- 
bard model with unit-matrix hopping studied here does 
not phase-separate in spite of having characteristics at 
n—1 very similar to those observed in the MC simula- 



tions of the JT-model at the same density. Having a 
FM/OAF state at n=l apparently is not sufficient for 
PS to occur .e2I This issue is conceptually important and 
will be addressed again in the next subsection when re- 
sults with a more realistic nondiagonal hopping matrix 
are analyzed. 

Due to the stability of the intermediate phases away 
from n=l it is not too surprising to observe a phase dia- 
gram at n=0.5 with similar characteristics to those found 
in Fig. 7(b). In Fig. 9(b), the unphysical region J>U' 
has only been analyzed briefly, just sufficiently to confirm 
that CO characteristics exist there, with relevant momen- 
tum fc=7r/2. In the regime J<U', the spin FM and in- 
commensurate phases are still stable, and T^(fc) presents 
a broad peak at fc=7r/2 compatible with the tendencies 
toward charge-ordering that appear at n—0.5. The prob- 
able pattern of orbitals here may have orbital "a" below 
"b" at site i, a mostly empty i-|-l site, the reverse orbital 
pattern at i+2, and another empty site at i-|-3, with this 
arrangement repeated in space. 



B. Nondiagonal Hopping Matrix 

The unit hopping matrix used in the previous subsec- 
tion is a simple and natural choice to gather qualitative 
information about the two-orbital problem. However, the 
studies of models designed for manganites need a more 
complicated hopping matrix, as shown in Sec. HI. In this 
subsection, the realistic hopping amplitudes along the y- 
direction are used to contrast the results against those 
obtained with the unit hopping matrix. The results are 
actually found to be dramatically different at densities 
different from unity. The organization of the subsection 
is similar as the previous one, i.e., the analysis starts with 
density n—1, establishing the main features of the phase 
diagram, and continues with ji^l with emphasis on the 
possible appearance of PS. 



1. Case ofn=l 

In Fig. 10(a), the spin structure factor S{k) is shown 
at [/'=30 for representative values of J. At J=10 and 
smaller, the signal was found to be clearly antiferromag- 
netic with a sharp peak at fc=7r. In the intermediate 
region, the spin structure is complex with some incom- 
mensurate characteristics, as exemplified by the result at 
J=12. For larger values of J, such as 15 and 20, the 
system becomes ferromagnetic or quasi-ferromagnetic. If 
J is increased further beyond the J—U' boundary, spin 
incommensurate structures have been observed, as in the 
study leading to Fig. 7(b) (results in this unphysical 
regime will not be discussed further here). In Fig. 10(b), 
the charge structure factor N{k) is shown at the same 
couplings used in Fig. 10(a). This quantity only devel- 
ops some nontrivial structure as the J—U' fine is ap- 
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proached. In this case indications of charge-ordering at 
k=TT are observed. In Fig. 11(a), the results correspond- 
ing to the orbital structure factor are shown. At small J, 
the order is uniform since T^(fc) peaks at k—Q. At J— 12, 
an incommensurate structure appears in the orbital sec- 
tor, similar to the results obtained for the unit-matrix 
hopping away from n—1, and related with the spin in- 
commensurability observed in S{k). For J =15 and 20, a 
robust orbital-staggered pattern is reached. 




FIG. 10. (a) Spin structure factor S{k) vs. k using the 
realistic hopping matrix along the a;-axis for a 20-site ID 
chain at n=l and (7' =30. The dashed, solid, dotted, and 
dotted-dashed lines correspond to J=10, 12, 15, and 20, re- 
spectively, (b) Charge structure factor N{k) vs. k for the 
same parameters as in (a). 

With the information obtained at C/'=30 and various 
J's, supplemented by results gathered at U'=5, 10 and 
20, a rough phase diagram can be constructed, shown in 
Fig. 11(b). In this case four regimes are identified. At 
small J, an AF/OF state appears which does not exist 
for the unit-matrix hopping. The reason is the following: 
For the realistic hopping matrix, orbital "a" has more 
mobility than the other. Then, to improve the kinetic 
energy the ground-state prefers to have those orbitals as 
the lowest-energy ones at every site. However, at n=l, a 
spin aligned orbital-uniform arrangement would not have 
any mobility due to the Pauli principle. Then, the op- 
timal situation is achieved with antiferromagnetic order 
in the spin sector. The presence of this state is the main 
difference at n=l between the results found for the unit 
hopping matrix shown in Fig. 7(b) and the results of 
Fig. 11(b). On the other hand, note that the AF/OF 
order does not take advantage of J since the use of the 
hoppings in this state leads to antiparallel spins at the 
same site and, thus, to an energy penalization propor- 
tional to J. For this reason, as the coupling J grows, a 
transition is expected to a state in close competition with 
the AF/OF one, namely, the FM/OAF state that has ap- 
peared several times in our investigations. In the latter 



J plays an important role since it favors spins alignment 
when two electrons visit the same site, leading to a spin 
ferromagnetic state. To optimize the kinetic energy, the 
optimal orbital pattern must be staggered, as discussed 
in the previous section when explaining Fig. 7(b). The 
interpolation between these two competing states is diffi- 
cult to predict and the computational work suggests that 
it proceeds through a complicated mixture of FM- and 
AF-orbital characteristics that lead to an overall spin and 
orbital incommensurate pattern. Only further studies in- 
corporating finite-size analysis will clarify if this regime 
is indeed stable in the bulk limit. 




FIG. 11. (a) Orbital structure factor T^(fc) vs. k for 
the same parameters as in Fig. 10(a). (b) Phase diagram for 
n=l for the realistic hopping matrix. The four phases are 
AF/OF, IC/OIC, FM/OAF, and CO, where "OF" denotes 
orbital-uniform and charge-disordered. The "CO" has not 
been characterized in detail besides its CO properties. 



2. Case ofn^l 

It is important to investigate how the phase diagram 
shown in Fig. 11(b) evolves as a function of density. To 
establish a connection with the results found in the case 
of the JT-phonons at large A, once again the FM/OAF 
state is here studied in detail since a state with simi- 
lar characteristics was indeed observed at n=l (Sec. Ill) 
in the mean-field studies with phonons, and the re- 
duction of the|-d|ensity led to a phase-separated regime 
in simulations .O Would the same occur in the purely 
Coulombic case with non-diagonal hopping? To gain in- 
sight, S{k) is shown in Fig. 12(a) at some representative 
densities using couplings U'=30 and J=15 that corre- 
spond to a point inside the FM/OAF phase of Fig. 11(b). 
It is interesting to observe an abrupt change in S{k) 
when n is varied from 1 to 0.9, with an incommensurate 
structure appearing in the latter. This incommensura- 
tion continues up to n=0.5. The charge structure factor 
N{k) also presents an abrupt change away from n=l, 
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with tendencies to charge ordering maximized at n=0.5 
(not shown). In addition, T^(k) drastically switches from 
an orbital-staggered pattern at n—l, to a uniform ar- 
rangement at nj^l as shown in Fig. 12(b). At n=0.5 and 
working with U'=10 and 30, it was observed that the 
spin-incommensurate and orbital-uniform characteristics 
are independent of J to a good approximation, as long 
as U'>J. 




FIG. 12. (a) Spin structure factor S{k) vs. k for the re- 
alistic hopping matrix, U'=30, and J=15 in the 20-site chain 
for several electron numbers. Incommensurate characteris- 
tics appear abruptly upon doping the n=l state, (b) Orbital 
structure factor T^(fc) vs. k for the same parameters and elec- 
tron number as in (a). Again an abrupt change from OAF to 
OF features is observed as n is reduced from unity. 



To understand the curious behavior reported in 
Figs. 12(a) and (b), the density vs. /x is plotted in 
Fig. 13(a) to search for the PS tendency. An anomalous 
behavior is observed at n=0.9, which has a tiny /x-window 
of stability, and thus, a large compressibility k (since k 
is proportional to dn/d/x). This behavior indicates a ten- 
dency toward PS in the data, and it may occur that small 
additions to the two-orbital Hubbard model such as JT- 
phonons or even the analysis of slightly larger clusters, 
may render the system truly phase-separated. To con- 
firm the presence of PS precursors, the local density rii 
is shown in Fig. 13(b) at n=0.9. Note that the DMRG 
method works with open boundary conditions, and rii is 
not necessarily equal to n at every site due to the lack of 
translational invariance. It is clear from this figure that 
large charge inhomogcneitics are present in the ground 
state, while at n=l (not shown) rii is almost uniform. 
Only the center of the chain has a density equal to the 
average one (0.9), while the ends of the chain have n close 
to unity, and at other sites near the center the density 
is smaller than 0.9. The minimum in the local density 
shown in Fig. 13(b) may correspond to a hole doped into 
the n=l system, that has developed polaronic character- 
istics. By symmetry, the other hole is on the other half of 



the chain. In Fig. 13(c), two holes generate two minima 
in the density (again, with the other two holes on the 
other chain half). These results show that tendencies to- 
ward charge inhomogeneities are present in this system, 
and precursors of PS are observed possibly in the form 
of polaronic behavior. The purely Coulombic model ap- 
pears to be at the verge of phase separation. 
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FIG. 13. (a) Density n vs. /j for the same parameter as 
in Fig. 12(a). Precursors of PS at density close to unity are 
observed, as discussed in the text, (b) Local density rii vs. 
site position i for 18 electrons on a 20 site chain. Only half 
the chain is shown, the results for the other half are simply 
obtained by reflection. The system is far from uniform in the 
charge sector, (c) Same as (b) but for 16 electrons. 

To search for more clear indications of phase- 
separation, the Coulombic interaction U' was reduced 
to 10, and J was fixed to 7. This still corresponds to a 
point at n=l inside the spin-FM orbital-staggered phase. 
The spin and orbital structure factors for several number 
of electrons were calculated in this case, and the results 
(not shown) have clear similarities with those obtained at 
U'=30. However, in this case now the analysis of the den- 
sity vs fi reveals strong phase separation characteristics 
between densities 1.0 and 0.7 (Fig. 14 (a)), qualitatively 
similar to those reported using JT-phonons. In Figs. 14 
(b) and (c), the local density away from n=l is shown. 
As in the case reported in Figs. 13 (b) and (c), strong 
oscillations reveal clear tendencies to phase-separation. 
Studies at U'=20 and J=12 but with the realistic hop- 
ping matrix have also been carried out as part of this 
effort. The results arc; very similar to those shown in 
Figs. 13 and 14. Summarizing, either clear phase separa- 
tion or a strong tendency to such phenomenon exists in 
the ID purely Coulombic model as long as the hopping 
matrix is non-diagonal. 
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C. Influence of Nearest-Neighbor Repulsion 



As discussed in Sec. IV, V should not be too large in 
the actual material, since if it were strong, the CS struc- 
ture would be destroyed. However, it is important to 
understand whether the results change or not with the 
inclusion of V using the same unbiased technique as in 
the previous subsection. 




FIG. 14. (a) Density n vs. for U'=W and J=7. Clear 
phase separation between densities 1.0 and 0.7 is observed, 
(b) Local density rii vs. site position i for 18 electrons on a 
20 site chain. Only half the chain is shown, the results for 
the other half are simply obtained by reflection. The system 
is far from uniform in the charge sector, (c) Same as (b) but 
for 16 electrons. 
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FIG. 15. (a) Charge structure factor N{k) vs. k at n—0.5 
with the realistic hopping for U' — IO and J— 6 in the 20-site 
chain. Two values of V as shown and the enhancement of 
charge-order is clearly observed, (b) The local density rii vs. 
i for the same parametre as in (a). Shown are results for half 
the lattice, the other half can be found by reflection. 



First let us consider density n=l, on a chain of 20 
sites and C/'=30. In this case, few modifications were 
observed compared with the results obtained for V^O. 
Namely, (i) the boundary between the spin-AFM and 
spin-incommensurate phases at small J, compared with 
[/', shifted toward a smaller J, and (ii) close to J—U' 
the charge correlations were enhanced substantially fa- 
voring the staggered pattern of charge between ni=0 and 
2, which is not penalized by V . Even at values such as 
J— 20, i.e. not too close to J=U', this enhancement was 
noticeable. However, since there are no indications of 
such a pattern in experiments for manganites at n=l, 
there is no need to analyze this region in more detail in 
the present study. 

At n=0.5, it is naively expected that the V^-term brings 
the CO state, although only precursors of CO behavior 
were observed at V=0 and this density as discussed in 
the previous subsection. In Fig. 15(a), N{k) is shown at 
U'=10 and J— 6 both at V=0 and 4. The region near 
k—TT is clearly enhanced by V as expected. The real- 
space density is shown in Fig. 15(b) and a clear charge- 
staggered pattern is visible. The spin structure factor 
(not shown) is still peaked at k=Tr/2 as at V=0, com- 
patible with a staggered spin-arrangement involving the 
occupied (even or odd) sites. The orbital structure factor 
still has a large uniform component, although it has also 
developed a broad low-intensity peak at k—n. Overall 
the characteristics of the state stabilized by V are orbital- 
uniform, charge-staggered with period of two lattice spac- 
ings, and spin-staggered over the occupied sites (period 
of four lattice spacings). Thus, only when a sufficiently 
large nearest-neighbor repulsion is included, the n=0.5 
CO-state similar to that found in the JT-phononic model 
is reproduced in the purely Coulombic model in one di- 
mension. This result seems in contradiction with the 
discussion around the CS structure observed in the half- 
doped material, since it was concluded that V is not the 
only origin of the CO-state in the manganite. However, 
as briefly discussed in Sec. HI C, if the zigzag ID chain, 
not the straight ID chain, is used in the calculation, the 
CO phase can be obtained even in the purely Coulombic 
model without V, although the (3x^ — r^)/(3?/^ — r^)-type 
00-phase cannot be obtained. From this viewpoint, a 
DMRG study carried out on the zigzag ID chain will be 
interesting, but this is left for the future. 

The analysis of the two-orbital model that contains 
only Coulomb interactions was instructive in several re- 
spects. For instance, the study has shown that a spin- 
ferromagnetic orbital-staggered (FM/OAF) phase ap- 
pears naturally in this context at n—1 for a variety of 
hopping amplitudes, result in excellent agreement with 
those observed for a purely JT model, also in agreemept 
with the previous ED studies of Coulombic models.tll 
Then, the FM/OAF phase is a robust feature of mod- 
els for manganites and approximations that attempt to 
describe these materials cannot neglect orbital-ordering. 
Another important result of this section has been the 
observation of PS tendencies with Coulombic-only inter- 
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actions. This was shown to occur in the form of robust 
precursors of this tendency in some regions of parame- 
ter space. However, to observe PS it seems necessary 
to use reahstic electronic hopping amphtudes in the case 
of Coulombic models. The presence of PS tendencies 
both using purely JT and Coulombic interactions con- 
firms the robustness of this feature, and it also shows 
that mixed-phase tendencies cannot be ignored in theo- 
retical studies of manganites. In addition, the similarities 
between models with JT-phonons or Coulombic repul- 
sions suggest that the technically much simpler studies 
that use only phonons to mediate the interaction between 
electrons are qualitatively correct, and likely capture the 
physics of more involved models where both interactions 
are included. 



ing the "Coulombic" and "JT-phononic" approaches to 
manganites as qualitatively different ways to carry out 
theoretical calculations is likely incorrect. 
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VII. DISCUSSION AND SUMMARY 

In this paper, the two-orbital model with the JT- 
phononic and/or the Coulombic interactions has been 
studied using a variety of techniques. Three points have 
been confirmed in this work, (i) The main properties of 
manganites can be reproduced successfully by the purely 
JT-phononic model even if the strong on-site correlation 
is not included explicitly, since the effect of such cor- 
relation can be renormalized into the effective electron- 
phonon coupling, (ii) In particular, in the mean-field 
level, the JT-phononic model can successfully reproduce 
the CE-type AFM phase with the charge-stacked struc- 
ture of the 3D cubic lattice at n=0.5, in excellent agree- 
ment with the MC simulations of Yunoki et al. Even 
if the nearest-neighbor repulsion V is introduced, this 
phase is not easily destroyed due to the key role played by 
the magnetic energy gain, regulated by Jaf, in the CE- 
vs. sCi?-type competition and the "topological" energy 
gain in the CE- vs. C-type competition. Particularly, it 
is stressed that the topology of the zigzag ID path is the 
key issue leading to the stabilization of the CO/00 state 
in the C-E-type structure. The on-site Coulombic model 
treated in the mean-field approximation and without JT 
phonons was also found to lead to charge- stacking due to 
the influence of Jaf. (hi) The purely Coulombic model 
behaves in many respects very similarly to the purely 
JT-phononic one and, in particular, it presents the phase 
separation tendency, especially when realistic hoppings 
are used. 

Summarizing, both approaches to the problem of man- 
ganites, based either on Coulomb repulsions or phonons, 
share common tendencies. This |Conclusion is in agree- 
ment with the recent observationllfl that the percolative 
character of transitions in manganites and its large mag- 
netoresistance effect arise from the competition between 
metallic and insulating phases in the presence of disor- 
der, independently on whether these phases are mainly 
generated by Coulombic or JT interactions. Our results 
have provided robust arguments suggesting that perceiv- 



APPENDIX A: RELATION WITH 
PHASE-SEPARATION THEORIES FOR 
CUPRATES 

Since a substantial portion of the paper is devoted to 
the issue of phase separation in models with Coulomb in- 
teractions, in this Appendix extra comments are provided 
regarding PS in models for high-temperature supercon- 
ductors (HTSC), issuGpwhich has been under discussion 
for almost a decadcPiPi. In particular, it is important 
to clarify the relation between HTSC phase separation 
and the PS phenomenon discussed here for manganites. 
In fact, some authors strongly believe that any doped 
correlated insulator should have PS, and in this respect 
the results for manganite models would be, a mere par- 
ticular case of a more general framework. E3 However, it 
has to be discussed in detail to what extend there is con- 
vincing theoretical evidence that indeed any doped insu- 
lator phase separates. If this is not clear, the relation 
between PS in cuprates and manganites weakens consid- 
erably. The discussion in this context also has to involve 
two other important aspects of the problem, namely the 
microscopic origin of PS, and its phenomenological con- 
sequences, particularly when extended Coulomb interac- 
tions are included. 

There are several differences between the PS phenom- 
ena proposed for manganites and cuprates. (i) The prop- 
erties of the two competing phases are not the same. In 
manganites, an undoped AFM-phase and a hole-doped 
FM-state are involved, while in cuprates it is an undoped 
AFM-state and a hole-doped paramagnet or supercon- 
ductor. In models for HTSC, ferromagnetism does not 
play an active role at realistic J/t couplings, while it is 
crucial in manganites. Nevertheless in Ref. ^ PS be- 
tween a G-type antiferromagnet and a ferromagnet was 
also discussed, within the framework of the very small 
J/t limit of the t-J model and assuming a fully saturated 
FM-state. The argument in Ref. |5^ can be applied di- 
rectly to the case of the one-orbital model for manganites 
to justify the presence of PS in this context. However, 
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note that the AFM-state that is doped in 3D real man- 
ganites has staggered spin order only in one direction, 
while it is ferromagnetic in the other two (A-type AFM 
order). In fact, PS can occur in two-orbitals models for 
manganites in ID and 2D without actually involving the 
AFM-state, but only having two competing FM-states, 
as discussed in the main tcxtO Here the orbital degree 
of freedom and JT phonons play the key role needed for 
PS, while they are not important in the cuprates. (ii) 
PS in models for manganites occur even in one-dimension 
in regions of couplings that are realistic, as described in 
this paper and previous ones, unlike t-J model results in 
the same dimension where PS only happens at large val- 
ues of J/t, and it does not occur at all in the one-band 
Hubbard model, (iii) The study of PS in models ana- 
lyzed in the infinite dimensional limit lead to different 
conclusions between cuprates and manganites, although 
this issue is still controversial. Fair is to say that the 
case of the D=oo Hubbard model is subtle and deserves 
a special discussion (D is dimensionality). In principle 
at D=oo there ace no indications of PS in the one-band 
Hubbard model.LJ However, it is not clear if this can be 
taken as a counterexample of PS in a Mott insulator, as 
recently remarked in Ref. ^ In most D=oo investiga- 
tions, the ratio J/t scales to zero like and as a 
consequence the physics of PS is difficult to address. In 
addition, when non-bipartite lattices are used, the AFM 
order is frustrated. Then, more work is needed in large 
dimension to clarify the effect of doping on correlated in- 
sulators. However, note that in calculations for the one- 
orbital Kondo model for manganites PS has been clearly 
observed even at D=oo (see Ref. 0). This suggest that 
important qualitative differences may exist between the 
PS phenomena in models for cuprates and manganites. 

Let us elaborate more on stripes in manganites and its 
topological vs. non-topological character. Whilp-stripe 
order has been claimed to exist in manganites J3 note 
that the regime of hole density in which this observation 
was made {x=2/3) is far from the low hole concentra- 
tion were the discussion of stripes in the cuprates occurs. 
These authors are not aware of experiments reporting 
stripes at small x in manganites, although charge in- 
homogeneities |Of various forms appear in several cases 
in this regime.EHl In addition, it is unc 
bi-stripe structure found by Mori et al. 
losical characteristics. Recent calculations by Hotta et 
aZ.E3 have shown that the manganite bi-stripes are better 
interpreted as arising from zigzag conducting chains run- 
ning perpendicular to the charge stripes once electron- 
phonon JT-couplings are switched-on. For this reason, it 
is premature to establish connections between stripes in 
cuprates and manganites. 

Regarding the issue of whether any doped correlated 
insulator produces PS, the following is our understand- 
ing of the current theoretical literature for the cuprates. 
The most clear manifestation of PS appears in the t-J 
model at large J/t, where the tendency to form pairs 
of holes to minimize the number of broken AFM bonds 
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is so intense that clusters of holes are formed instead 
of individual pairs. Here the attractive potential energy 
among hole carriers originating in the AFM background 
dominates over the kinetic energy, that tends to spread 
particles apart. As J/t is reduced it has been a matter 
of much controversy whether the PS effect survives in 
the realistic small J/i—regime of the t- Jj-nmdel. While 
ED of small clusters,EEl|-MC simulations j£3'LiJ and high- 
temperature expansionsE3 suggested that at J/t oLiprder 
unity the effect would disappjear, other argumenta^il and 
further computationalworkE^I opened the possibility for 
PS to jfixist at all values of J/t. Very recent DMRG 
studiesEj for ladders of incw;asing number of legs show 
that previous calculationa£3 may not have been suffi- 
ciently accurate and the new results suggest that PS in- 
deed only occurs at intermediate and large values of J/t 
in the 2D t-J model. This same conclusion becomes more 
clear once extra hopping amplitudes are added to the 
model. In this case the substantial hole mobility induced 
by the extra hoppings shifts the PS regime_tp values of 
J/t larger than in the pure t-J model case.Lj Then, the 
proposal that any correlated insulator when hole doped 
should become phase separated is still not confirmed us- 
ing unbiased techniques in simple models for cuprates. 

Regarding the phenomenological aspects of the PS 
regime, the issue of microscopic PS that may appear 
in manganites has certainly been—discussed before by 
Emery and Kivelson for cuprates O In this context the 
Coulomb interactions break into small pieces the macro- 
scopic clusters of the two phases in competition, since 
they have different electronic densities. Stripe patterns 
emerged froro-calculations carried out mainly close to the 
atomic limit j£j where the attraction leading to PS plus 
the Coulomb repulsion are in competition. In the con- 
text of Nudear Physics similar patterns have also been 
discussed.Ell In addition, Nagaev studied the formation 
of finite size clusters of one phase embedded into tha 
other, mainly for antiferromagnetic semiconductors .Ea 
Then, the simple picture of a stable state formed by 
small clusters of the competing phases, somewhat sim- 
ilar to the CDW pattern obtained in the ID calculations 
with nearest-neighbor repulsions, is certainly common to 
manganites and cuprates, and it has been described—in 
the context of the Frustrated PS scenario for HTSCo 

In short, the discussion presented in this Appendix sug- 
gests that the PS phenomena in models for manganites 
and cuprates, while sharing the general common ingre- 
dients of any PS regime, are different in origin and they 
must be considered separately in their study. 
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